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Abstract We present Binomials, a package for the computer algebra system Macaulay 2, 
which specializes well known algorithms to binomial ideals. These come up frequently 
in algebraic statistics and commutative algebra, and it is shown that significant speedup 
of computations like primary de composition is possible. While central parts of the im- 
plemented algorithms go back to lEisenbud and Sturmfels (Il996l) . we also discuss a new 
algorithm for computing the minimal primes of a binomial ideal. All decompositions make 
significant use of combinatorial structure found in binomial ideals, and to demonstrate the 
power of this approach we show how Bino mials was used to c ompute primary decomposi- 
tions of commuting birth and death ideals of lEvans et a yielding a counterexample 
for a conjecture therein. 

Keywords algebraic statistics • binomial ideals • commuting birth and death ideals • 
computational commutative algebra • primary decomposition 



1 Introduction 



A monomial ideal is an ideal generated by monomials, a binomial ideal is one whose gen- 
erators can be chosen as binomials. A pure difference ideal is an ideal whose generators are 
all differences of monic monomials. For monomial ideals, central concepts Uke Grobner 
bases, irreducible and primary decompositions, etc. can be defined directly on the expo- 
nent vectors of the monomials generating the ideal. In this sense the whole theory is very 
combinatorial. For binomial ideals the si tuation is more complicated, bu t essentially it can 
be made combinatorial too. Starting with lEisenbud and Sturmfelsl ll 19961) the combinatorial 
theory of binomial ideals has developed into a branch of combinator ial commutative alge- 
bra w hich has many connections to different areas of mathematics dMiller and Sturmfels 
2005h . 

The interest in binomial ideals is motivated by the frequency with which one encounters 
them. For instance, commutative semigro up rings are e xactly the quotients of polynomial 
rings by pure difference binomial ideals tailmeij ll984l) . Toric id eals, which are binomial 
prime ideals, are the defining ideals of toric varieties as defined bv lFultonl(ll993h . This fact 
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is central in the field of algebraic statistics, where closures of discrete exponential fami- 
lies, such as graphical or hierarchical mo dels, have been recognized to be nonnegative real 
parts of toric varieties jGeiger et a l l2006h . Also binomial ideals which are not prime occur 
there. Conditional independence models are defined through a set of polynomial equations 
in the elementary probabilit ies, and studying primary decompositions of the corresponding 
id eals is of natural interest (Drton e^a il l2009l: iFinkl l I2OO9I: lHerzog"etall2009l) . For instance, 
as lEisenbud and Sturmfelsi (il996i) have shown, the minimal primes of binomial ideals are 



essentially toric ideals, and therefore a conditional independence model is a union of expo- 
nential families. Knowing the primary decomposition, a piecewise parameterization of the 
model is instantly available. 

This paper deals with the polynomial ring ]k;[;ici, . . . ^x„] over a field k of characteristic 
zero. Choices for k are the rationals Q, their cyclotomic extensions Q(^/), or the complex 
numbers C. Primary decompositions of binomial ideals are not necessarily binomial as is 



easily seen on the ideal (x — 1), which over Q decomposes as {x — 1) H (x 



1). If k 



is algebraically closed, however, binomial primary decompositions exist. When speaking 
of primary decompositions in this paper, we always mean primary decomposition into bi- 
nomial ideals, and we have to extend the coefficient field where needed. For the software 
package we have restricted even further: We consider only pure difference binomial ideals. 
In that case, the primary decompositions into binomials will be shown to exist with coef- 
ficients in cyclotomic extensions of Q. In m any applications it suff ices to study this case. 
Examples include the semi-graphoid ideal (iHemmecke et a conditional indepen- 

dence ideals, commuting birth and death ideals of Section|2l and almost any other binomial 
ideal considered in algebraic statistics. 

This paper is structured as follows: In Section [TTI we study a systematic way of ap- 
proximating binomial ideals by cellular binomial ideals. Then in Section [L2l we give an 
algorithm for finding the solutions of zero-dimensional pure difference binomial ideals and 
apply it to saturation of partial characters. In Section [T3] we give a new algorithm for com- 
puting the minimal primes of a binomial ideal. Section|2]contains results on large primary 
decompositions that have been carried o ut with our softwar e Binomials. We show a coun- 
terexample to Conjectures 5.3 and 5.9 in Evans et al ( 2010l) . Finally, Section [3] concludes 
the paper with future research directions. 



Th roughout the paper we use notation that tries to coincide with that of lEisenbud and Sturmfels 
(Il996h . We assume familiarity with basic notion s of commutative algebr a. A very pedagog- 
ical introduction is the book of Cox et al (Il996l) . while .Eisenbudid 1995b covers everything 
from the very basics to current research topics. In keeping with the introductory nature of 
this work, each of the following sections contains examples of how to do the discussed com- 
putations with the help of Binomials. These examples are thought of as a motivation and 
do not cover all of the functionality that is implemented. They are produced with version 
0.5.4 of Binomials. The reader is encouraged to download the package, use it, and report 
experiences to the author An online help is integrated. 

Example 1 (Installation) Binomials and an auxiliary package for cyclotomic fields, called 
Cyclotomic, are available under the URL: 

|http : / / personal- homepages . mis . mpg . de/kahle/bpd/| (1) 

It is recommended to install the latest version of Macaulay 2 ( Eisenbud et all 2001 ) before 
using Binomials. To get started, run Macaulay 2, then load the package with 



il 



load 



Binomials . m2 " 



The additional packages FomTiTwo and Cyclotomic are needed. The first is included in 
Macaulay 2 as of version 1 .2, while the latter can be obtained together with Binomials. To 
make the documentation available the package should be installed: 

12 : InstallPackage ("Binomials", RemakeAllDocumentatlon=>true) 

After running this, help can be accessed with 
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i3 : help "Binomials" 



1.1 Cell Decompositions of Binomial Varieties 

Our analysis of a binomial variety starts with the decomposition of k" into the 2" algebraic 
tori interior to the coordinate planes. Each of the coordinate planes is defined by a subset 
^ C {1, . . . ,n} of the indeterminate's indices. We denote the algebraic torus corresponding 

ioShy 

(k*)'^ := {(xi, . . . ,x„) e k" : Xi 7^ 0,/ e and x^- = 0,Vy ^ <f } . (2) 

Geometrically, for a binomial ideal / C k[xi , . . . ,x„], we study cellular decompositions. 
Their components are the intersections of primary components which have generic points 
in a given cell (k*)*^. The central definition is 

Definition 1 A proper binomial ideal / C h[xi ,x„] is called cellular if each variable x, 
is either a nonzerodivisor or nilpotent modulo /. 

In this paper a variable is always a variable in a polynomial ring, random variables are 
not mentioned explicitly. Primary ideals / are cellular as every element of k[xi , . . . ,x„]// is 
either nilpotent or a nonzerodivisor. The following explicit representation of cellular ideals 
is only a reformulation of the definition but useful in many ways. 

Lemma 1 A binomial ideal I C k[xi,...,x„] is cellular if and only if there exists a set 
^ C { 1 , . . . , «} of variable indices of k[xi , . . . ,x„] such that 

2. For every i ^ S, there exists a nonnegative integer di such that the ideal (xf : i ^ S 
is contained in I. 

We call the set S' the cell indices and the variables {x, : i € S'}, which are exactly the nonze- 
rodivisors modulo /, the cell variables. We denote by M{S') the ideal generated by the non- 
cell variables, i.e. the variables {x, : / ^ For any vector d — {di)i^g of natural numbers 

we denote M{SY '■= (x^l' '-i ^S^. With this notation, another useful representation of cel- 
lular ideals is given by the following Lemma. In lEisenbud and Sturmfels (1 19961) the ideal 



on the right hand side of (O is denoted 7^''' . 

Lemma 2 A binomial ideal I is cellular if and only if there exist a set <S (- {1 , . . . , n} and 
an exponent vector d, such that 

I=(^[l + M{£Y'):(^llx}j y (3) 

Radicals of cellular binomial ideals have a nice combinatorial structure, defined by the 
set (o, and a partial character, which we introduce next. For this let ^ C {1, ...,«} be 
any nonempty subset of the indices of variables and define the shorthand k[(o'] :— k[x, : / 6 

Definition 2 A partial character is a pair a), consisting of an integer lattice C 
Z"^ and a map CJ : ^ — > k*, that is a homomorphism from the additive group ^ to the 
multiplicative group k*. For each integer lattice ^ C Z*^, we define its saturation 

Sat(Jf ) |m e if : dm e ^ for some t/ e Z | . (4) 

A lattice C TL^' is called saturated if it satisfies J£ ~ Sat(^). A partial character cj) 
is called saturated if ^ ~ Sat(^), and it is called a saturation of a partial character 
{^', a'), provided that ^ = Sat(^0 and C7'(/) = (7(/), V/ e 
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Often it is convenient to denote by L an integer matrix having the lattice ^ as its right 
image ^ := {Lm :mElf^. Thus, the columns of L span the lattice, and we abuse notation 
speaking of the partial character (L, a) in this case. To each partial character a) we 
associate a lattice ideal: 

I^{a) := (x"''' - a{m)x"'' : m e C (5) 

Here we have decomposed m G l/ into its positive part mf := max{OT;,0}, and negative 
part , so that m = ;«+ — m^. We also have used monomial notation x™ := Ylieg^T ■ ^'^ 
the notation of (|5]l, the lattice is alwa ys implicitly understood from a . 



It follows from Theorem 2.1 of lEisenbud and Sturmfelsl (Il996h that a lattice ideal is 



prime if and only if its partial character is saturated. More generally, all associated primes 
of a lattice ideal arise from saturations of its partial character A nice characterization is that 
a proper binomial ideal / C ]k;[xi , . . . ,x„ ] is a lattice ideal if and only if / = (/ : (ITLi -*^0°°)- 
This fact can be used to compute a minimal generating set of a l attice ideal when only 



the pa r tial character is giy e n, a problem considered for in stance in Ho^ten and Sturmfelsl 
(ll995b : lBigatti et all (fl999l) : lHemmecke and MaUdnl ( |2009|) . 



A cellular binomial ideal is a lattice ideal on a subset of the variables. For instance, 
it follows from Lemma|2]that radical cellular binomial ideals / C k[xi , . . . ,jc„] are of the 
form / = M{S') +I+{o) for some partial character (L, o) on l/ . Now, assuming that k is 
algebraically closed, the associated primes of M{S') +I+{a) are given by 

P, = M(#)+/+(t), (6) 

where T runs through all saturations of a. In particular a radical cellular binomial ideal 
is equidimensional. If k is not algebraically closed, it may contain only some, or even 
no saturations of {L,a). In Section [T31 we give an algorithm that computes the minimal 
primes of a binomial ideal by directly finding a cellular decomposition of the radical of / 
into radical cellular ideals. 

If the monomials in a cellular binomial ideal / are of higher order, then we only have 
that /n k[<#] is a lattice ideal. However, the associated primes might have partial characters 
supported on different lattices. The key theorem for computing associated primes of cellular 
binomial ideals is 



Theorem 1 dEisenbud and SturmfeS (Il996l) . Theorem 8.1) Let I C k[jci , . . . ,Xn] be a 



cellular binomial ideal on the cell variables S. Let P = M{S') +/+((7) be an associated 
prime ofL then there exists a monomial x"^ in the variables not in S and a partial character 
T on whose saturation is (7, such that 

(/:x"')nk[^]=/+(T). (7) 

Note that the associated primes of a cellular binomial ideal are cellular binomial ideals 
for the same cell variables. To compute them, one considers all quotients of / modulo 
the standard monomials in the variables outside S". There are only finitely many, as / is 
cellular and contains M{(^y' for some nonnegative integer vector d. This theorem reduces 
the computation of associated primes to cellular decomposition and saturation of partial 

characters. 

We now review an algorithm for computing cellular decompositions due to Oieda and Sanche3 



(I2OOOI) . [t is based on the following approximation scheme for arbitrary ideals in any 



Noetherian ring: 

Lemma 3 dEisenbud and Sturmfelsl ( Il996l) . Proposition 7.2) Let I be an ideal in a Noethe- 



rian ring S and g € S such that (I : g) = (I : g°°). Then 

L I^{I:g)n{I+{g)). 

2. Ass(5/(/ : g))nAss{S/{I+ {g))) = 0. 
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3. A minimal primary decomposition of I consists of the primary components of {I : g) and 
those primary components of I + (g) that correspond to associated primes of I. 

Given any noncellular binomial ideal /, we can find a variable Xi that is a zerodivisor but not 
nilpotent modulo /. A power 5 > of that variable satisfies the conditions on g in Lemma 
|3]and we can write 

/=(/:4)n(/+(4)), (8) 
where the ideals on the right hand side are both binomial and properly conta ining/. This can 



be turn ed into a simple algorithm for cellular decomposition, formulated bv lOieda and Sanchez 



(l2000h . The authors also provided an implementation in Macaulay 2, parts of which are 



used in the Binomials package. 

Algorithm 1 (Cellular Decomposition) Input: /, a binomial ideal. 
Output: A cellular decomposition of /. 

1 . If / is cellular, return /. 

2. Choose a variable that is a zerodivisor but not nilpotent modulo /. 

3. Determine the power s such that (/ : x^- ) = (/ : 

4. Iterate with (/ : jcj) and /+ (x-). 

Step [T] is carried out as follows. First determine the nilpotent variables by checking for 
which X; one has (/ : x°°) = k[xi , . . . ,x„]. Denoting the remaining variables' indices as S', I is 
cellular iff (/ : (IIig^^;)") = ^- Termination of Algorithm[T]is ensured since ]k;[xi , . . . ,x„] is 
Noetherian and the two ideals (/ : x^ ) and /+ (xj) properly contain /. Correctness follows 
from Lemma [3] Note also that cellular components of pure difference binomial ideals are 
pure difference binomial ideals. 

Example 2 ( Cellular Decomposition ) We study an ideal from lEisenbud and Sturmfelsl(ll996h . 



Let 5 = Q[xi , . . . jXg] and / = (xix| — X2X5,XjX3 — xfx2,X2X^ — X3X5). 

11 : S = QQ[xl,x2,x3,x4,x5] ; 

12 : I = ideal (xl*x4''2-x2*x5'2, xl " 3*x3 " 3-x4 " 2 *x2 " 4 , x2 *x4 " 8-x3 " 3*x5 " 6 ) ; 

13 : toStrlng BCD I 

03 = { Ideal (xl *x4 " 2-x2 *x5 " 2 , xl " 3*x3 " 3-x2 " 4 *x4 " 2 , 

x2"3*x4''4-xl'2*x3'3*x5"2, x2 " 2 *x4 " 6-xl*x3 " 3*x5 " 4 , 
x2*x4"8-x3"3*x5*6) , 
Ideal (xl"2,xl*x4"2-x2*x5'2, x2 " 5 , x5 " 6 , x2 " 4 *x4 ' 2 , x4 " 8 ) } 

14 : ap - blnomlalAssoclatedPrlmes I; toStrlng ap 

04 = (Ideal (xl*x4''2-x2*x5'2, xl " 3*x3 " 3-x2 " 4 *x4 " 2 , 

x2"3*x4"4-xl"2*x3"3*x5''2, x2 " 2 *x4 " 6-xl*x3 " 3 *x5 " 4 , 
x2*x4"8-x3''3*x5'6) , 
Ideal (x2,x5,x4,xl) } 

15 : Intersect (ap#0,ap#l) == I 

05 = false 

15 : blnomlalRadlcal I == Intersect (ap#0,ap#l) 

06 = true 

17 : IsCellular (ap#0, returnCellVars=>true) 

07 = {xl, x2 , x3, x4, x5 } 

18 : IsCellular (ap#l, returnCellVars=>true) 

08 = {x3} 

In this listing we have suppressed some output. First we compute a cellular decomposition 
with BCD. It has two components. The first ideal is the toric ideal (/ : (n"=i-*^i)°°)' which is 
prime. It is a general feature of the implementation of Binomials that, when the input has 
no monomial generators, the first ideal of the output of cellular and primary decomposi- 
tions, as well as minimal and associated primes, is always the toric ideal. We also compute 
the associated primes. The second one is embedded, and we confirm that / is not radical. 
Note also that the binomial generator X1X4 — X2X5 in the second cellular component reduces 
to zero as soon as one takes the radicals of the monomials. Finally we confirm that the 
associated primes are cellular, and show the set of variables with respect to which they 
are cellular, using isCeiiuiar with the option returnCeiivars. The cell variables could 
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have been computed directly together with the cellular decomposition by running the long 

version binomialCellularOecomposition, again with the option returnCellVars Set tO 
true. 

Theorem[T]shows that saturation of partial characters is a crucial ingredient for comput- 
ing associated primes of a binomial ideal. We therefore study the properties of saturations 
of partial characters. In the current implementation of Binomials any operation that needs 
extension of the coefficient field of the polynomial ring is only implemented for pure differ- 
ence binomial ideals. It will be shown that in this case cyclotomic field extensions suffice. 



1.2 Solving Pure Difference Binomial Ideals 

In this section we give a fast algorithm for solving pure difference binomial ideals of di- 
mension zero. It is not surprising that such a procedure utilizes only the exponents of the 
generators. We denote by ^/ the primitive l-th root of unity exp { ^ } G C. The field exten- 
sion of Q that is obtained by adjoining such a root of unity is called a cyclotomic field and 
denoted by Q(i^/). It can be constructed by taking the quotient of a univariate polynomial 
ring modulo the principal idea l generated by the minimal polynomial of the cyclotomic 
polynomial (lHungerfordll974 , Chapter V). 



Proposition 1 Given a zero-dimensional pure difference binomial ideal I, there exists a 
primitive root of unity ^/ such that all complex solutions of I are contained in the cyclotomic 
field Qi^,). 

The proof is given after the following Lemma, also of interest for the implementation. 
Lemma 4 The complex solutions of the univariate equation 

^ = C (9) 

are given by the following roots of unity 

■^0 — bmn' — hmn ' • • • I — S'nn ■ K'^'JJ 

Proof The xq, . . . ,jc„_i are n distinct roots of (|9]l, which is of degree n. □ 

Proof (Proposition\J} The standard method of reducing a multivariate problem to a uni- 
va riate prob l em ap plies. The general framework is described for instance in Chapter 3 
of lCox et all (Il996l) . Choose an elimination term order, such as lexicographic order, and 



compute a Grobner basis of /. This Grobner basis consists of pure difference binomials 
since all S'-polynomials are pure difference binomials. Furthermore, at least one of the bi- 
nomials of this Grobner basis is univariate as / is zero-dimensional and we have chosen an 
elimination order. The solutions of this univariate equation exist in a cyclotomic field by 
Lemma |4] We continue to extend the partial solution that we have found, substituting the 
variable for its value in the remaining elements of the Grobner basis. We obtain a univariate 
equation in another variable. The final solution exists in the cyclotomic field containing all 
the roots of unity that are encountered in the course of the algorithm. □ 

Of course, the procedure that was just described is also valid for other fields k. In the 
general case, field extensions have to be carried out by computing the minimal polynomial 
of the element to be adjoined and one has to do computations over the algebraic numbers. 
While possible in principle, this quickly becomes infeasible in practice, since both the 
computations become lengthy and it becomes more and more tedious to produce output in 
a human-readable form. 

We are now ready to formulate the algorithm for computing the variety of a zero- 
dimensional pure difference binomial ideal. The first thing that needs to be accounted 
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for is the possibility of as a solution, potentially with multiplicities. We take care of 
this by means of cellular decomposition. Each cellular binomial ideal / can be written as 
/ = ((/ + M((S')'') : (n,G^-^()°°)' and /nk[(?] is a lattice ideal. The solutions of / take the 
value zero at the variables outside S and each solution has a multiplicity of H;^^ di, where 
the di need to be chosen minimal. 

Algorithm 2 (Solving pure difference binomial ideals) 

Input: A zero-dimensional pure difference binomial ideal /. 

Outputs: The root of unity that needs to be adjoined to Q and the list of the solutions of /. 

1 . Compute a cellular decomposition of /. 

2. For each cellular component: 

(a) Set the noncell variables to zero, and determine the product D = Ylit^sdi of the 
minimal powers of the noncell variables. 

(b) Compute a lexicographic Grobner basis and solve the lattice ideal of the cellular 
component, adjoining roots of unity where necessary. 

(c) Save each solution D times. 

3. Compute the least common multiple m of the powers of the adjoined roots of unity and 
construct the cyclotomic field Q{^,„). 

4. Output the list of collected solutions as elements of Q(^,„). 

This algorithm is the main ingredient for saturating partial characters, which we treat after 
an example. 

Example 3 (Solving Pure Difference Binomial Ideals) We solve a simple pure difference 
binomial ideal to introduce the syntax. 

11 : S = QQ[x,y,z]; 

12 : I = ideal (x" 2-y , y " 3-z , x*y-z ) ; 

13 : binomialSolve I 

BinomialSolve created a cyclotomic field of order 3 . 

03 = {{1, 1, 1} r {- ww_3 - If ww_3, 1}, {ww_3, - ww_3 - 1, 1}, 

{0, 0, 0}, {0, 0, 0}, (0, 0, 0}} 

14 : degree I 

04 = 6 

In the implementation generic names consisting of ww and the order are assigned to roots of 
unity. Note that the square of the third root of unity ww_3 is represented as -ww_3-i by means 
of its minimal polynomial over Q. A cellular decomposition reveals that this ideal has 
two components, one of which is of degree 3 with associated prime {x,y,z)- The function 
binomialSolve outputs the Solutions with the correct multiplicities. If this is not desired, 
one can pass to the radical before solving, or directly compute the binomial minimal primes. 

Saturations of partial characters exist only over algebraically closed fields. This is evi- 
dent for instance from the partial character ((2), 2 i-t- —1), consisting of the rank 1 lattice 
spanned by the integer 2, and the character that maps 2 to — 1 G C. The saturations are pairs 
(Z, t), that satisfy t(2) = t(1)^ = —1. This example is merely a combinatorial version of 
factorizing the polynomial + 1, which is the same as performing the primary decompo- 
sition of its principal ideal. The following algorithm to saturate a partial character is the 
general version of the example's principle. 

Algorithm 3 (Saturation of a partial character) 

Input: A partial character (L, a), where L is a matrix whose columns are minimal genera- 
tors of a lattice in l/. 

Output: All distinct saturations (Sat(L), T,), ; = I,...,«. 
I. Compute the saturation L' := Sat(L). 
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2. Express the generators of L in terms of the generators of L' , by solving the matrix 
system 

L^L'K, (11) 

for the square matrix K = {kij)ij=\,...,r, where r := rk(L) = rk(L') denotes the rank of 
the lattices. 

3. Write I'j, and kj for the columns of L, L', and K, respectively. Introduce new vari- 
ables T, := = 1, . . . ,r, for the values that T takes on the columns of L'. Using 
again monomial notation t" : = 11!= 1 1/' compute the following zero-dimensional lat- 
tice ideal in Q[ti , . . . , Tr] 

J ■■= (^{t"' cy{Ij)T'J : y = 1, ...,/-): (Ot,) , (12) 

for the given values cj(/j). 

4. Solve J (over a suitable extension of Q) and output L' together with the list of solutions 
ofy. 

Proof (Correctness) Computing the saturation of a lattice should be viewed as an inte- 
ger valued analogue of taking the orthogonal complement twice. It can be carried out in 
Macaulay 2, for instance by computing the minimal syzygies of the syzygies among the 
generators of L. The coefficient matrix K that solves the system (fTTI) exists and is unique 
over Z, as L is a sublattice of L' and we assumed that the columns of L' are a minimal set of 
generators of the corresponding lattice. The ranks of L and L' coincide by definition. The 
ideal J is constructed as follows: For each generator Z of L we get a relation / — L' ■ k, to 
which we apply the homomorphism T, remembering that T and (7 are required to coincide 
on the generators of L. The entries of K are integers, thus we get the Laurent binomial ideal 

(oilj)~YlT^-^:j=l,...,ry (13) 

whose intersection with Q[ti , . . . , T,-] is exactly /. That J is zero-dimensional follows since 
the qu otient L'/L is a finite group. For details see Corollary 2.2 in .Eisenbud and Sturmfelsl 
(Il996h . □ 



The number of distinct saturations equals the order of the finite group Sat(L) /L, which 
can be computed by diagonalizing the matrix L, representing the inclusion Z'' — !■ Z^'. The 
Macaulay 2 command for this diagonalization is smithNormaiForm. Finally, for computing 
primary decompositions of pure difference binomial ideals we only need to solve such 
ideals during the saturation. 

Proposition 2 The saturation of a partial character that occurs during primary decompo- 
sition of a pure difference binomial ideal involves only solving pure difference binomial 
ideals. 

Proof Any cellular component of a pure difference binomial ideal is pure difference again. 
So we can assume that / is cellular Now, each partial character consists of a lattice and 
the constant map / 1— > 1. Therefore the ideal J in Algorithm [3] is a pure difference binomial 
ideal. □ 



1.3 Minimal Primes of Binomial Ideals 

In this section we describe a new algorithm for computing the minimal primes of a binomial 
ideal. It is based on a variant of cellular decomposition, given in Algorithm[T] As we have 
seen previously, the associated primes and thereby the minimal primes of a binomial ideal 
come in groups, associated to the cellular components of /. Our approach is to directly 
compute a cellular decomposition of the radical of /. 
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Algorithm 4 (Minimal primes of a binomial ideal) 

Input: A binomial ideal / C h[xi ,x„] . 
Output: The binomial minimal primes of /. 

1. Determine whether / is cellular. 

(a) If yes, compute the radical (/n k[(f]) +M((?) =M{S') +/+(c7) and its partial char- 
acter (L, cr). Compute the saturations of CJ and save the ideals 

P(t,-)=M(^?)+/+(t,-). (14) 

(b) If not, determine a variable x, that is a zerodivisor, but not nilpotent modulo /, and 
iterate with the ideals /+ (x,) and (/ : x°°). 

2. From all primes collected, remove redundant ones to find a minimal prime decomposi- 
tion of Rad(/). 

Proof (Termination and correctness) Termination of this algorithm follows as the ambi- 
ent ring is Noetherian and / + (x,) and (7 : x°°) strictly contain I . The radical of / is the 



intersection of the ideals in (4.2) of lEisenbud and Sturmfelsl (Il996h . We encounter a 



decomposition of Rad(/) into such ideals in the course of the algorithm, as the iteration 
is ultimately producing cellular components of the r adical of /. Thus, like in their Algo- 
rithm 9.2, correctness has been proved in Section 4 of lEisenbud and Sturmfels (Il996l) . For 



cellular ideals the minimal primes have the form (fl4] i. and the collection of all minimal 
primes of all cellular ideals contains the minimal primes of the original ideal by Lemma[3] 

□ 



This algorithm differs from the cellular decomposition algorithm only in the recursion step, 
where we continue with / + (x,) instead of / + {xf). In this way we do not achieve a decom- 
position of /, but only of the radical of /. Fortunately, this algorithm can be significantly 
faster than cellular decomposition since adding variables, instead of higher powers of vari- 
ables, allows the Grobner basis engine to do more simplifications during the computation. 

Example 4 (Binomial Minimal Primes) We continue where we left off in Example|2] 

il6 : toString binomialMinimalPrimes I 

ol6 = {ideal (xl*x4"2-x2*x5"2, xl " 3*x3 " 3-x2 " 4 *x4 ~ 2 , 

x2"3*x4~4-xl"2*x3"3*x5~2, 2 " 2 *x4 ~ 6-xl *x3 " 3 *x5 ~ 4 , 

x2*x4~8-x3~3*x5"6) } 

The result consists only of the toric ideal, confirming that the monomial prime is embedded. 
Although not visible from the output, the second associated prime was not computed on 
the way to this result. In particular, the minimal primes are not extracted from a list of 
associated primes. 



1 .4 Primary Decomposition 



The ori ginal primary decornpositio n algorithm of lEisenbud and Sturmfelsl (Il996h was re- 
fined bv lOieda and Sanchez (I2OOOI) . The computation starts with a cellular decomposition, 
a first approximation of primary decomposition. It is interesting to identify cases in which 
the cellular decomposition is already a primary decomp osition. Results in this direction are 
contained in Eisenbud and Sturmfels' paper, and also in lAltmannl(l2000l) . Note that in these 
cases a pure difference binomial ideal has a primary decomposition into pure difference 
binomial ideals, which is in particular independent of the coefficient field. 

In the general case, for each cellular component the associated primes need to be deter- 
mined. Then finding the primary component can be achieved as follows. From an associated 
prime P of a cellular binomial ideal /, extract the "binomial part" P'''^ — Pr\h.[S]. Then 
I + pW has P as its unique minimal prime. Computing the primary component over P is 
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carried out by means of a localization operatio n called Hull, remov i ng the embedded pri- 
mary components of I + P^''\ The refinement of Oieda and Sanchez ( 2000l) is to show that 
/ + f suffices in this procedure, while Eisenbud and Sturmfels originally suggested to 
add a sufficiently high mono mial power A combinato rial description of the resulting pri- 
mary components is given in 
these results for computation 



Dickenstein et all ( 2008 ). however, it seems difficult to use 



A few remarks on primary decom positions in Eisenbud and Sturmfels ( 19961) and lOieda and Sanchez 

(l2000l) are necessary. Corollary 6.5 of lEisenbud and Sturmfelsl ( 119961) shows that Hull(/) is 
a binomial ideal if / is a cellular binomial ideal. This corollary is used in the proof of Theo- 
rem 7.1' to deduce that Hull(/?,) is binomial, where /?, is the sum of a monomial ideal and 
/ + P*^*' from above. However, it is not checked whether /?, is in fact cellular, as required 
by the corollary. Example |5] shows a noncellular Rj that arises in the decomposition of the 
ideal of adjacent (2 x 2)-minors of a generic (5 x 5)-matrix. The computations necessary 
to check the example can be carried out easily with Binomials. 

Example 5 In the ring Q[a,b, . . . ,o] consider the ideal 

/ = (/« — ko,lt7i — jo, km — jn,l^,kl,jl,k^,jk,ik — hi, 
fk-cl, f , - gl , hj -gk,fj-al,cj-ak, fh - ci, 
fg — ai,cg — ah,f^,cf,af,ce — bf,ae — df,c^,ac,ab — cd,a^) . 

This ideal is cellular with respect to S' = {b,d,e,g,h,i,m,n,o}, and has four associated 
primes, which are pure difference. The binomial part of the unique minimal associated 
prime is 

P^''-* — {in — ho , im — go , hm — gn). 

Then I+P'-''^ has two cellular components whose sets of cell variables are S' and {b,d,e,m,n,o}, 
respectively. 



Using Theorem 7.1', in Algorithm 9.7 of lEisenbud and Sturmfelsl (1 1996h it is asked to com- 
pute Hull(P,), using Algorithm 9.6. This however, requires a cellular ideal as its input. The 
algorithm can be corrected easily since the operation Hull is called only for ideals whose 
radical is prime. The associated primes of such an ideal have the radical as their unique 
minimal element, and as Hull removes embedded primary components, instead of Hull(/?,) 
we can compute Hull(gi) of any other ideal Qi D Ri that has the same minimal prime. In 
particular we can c hoose Qj = (Rj : iY\ifz,ffXi)°°), th e "cellularization" of Rj. Summarizing, 



in Algorithm 9.7 of lEisenbud and Sturmfels ill 9961) Step 3.3 should be replaced by 



3.3' Compute Hull (/?,■ : (11(6^?-*^')°°) using Algorithm 9.6. 

Unfortunately, also i n The orem 3.2 of lOieda and Sanchez! (|2000|) . Corollary 6.5 of 



Eisenbud and Sturmfelsl (1 19961) is used to deduce that Hull(/+ (f n ]k;[(l'])) is binomial and 
primary. Again, this is wrong as /+ (f n ]k;[<o']) is not necessarily cellular The result can be 
saved by first cellularizing as explained above. The implementation in Binomials incorpo- 
rates these modifications and is demonstrated next. 

Example 6 (Binomial Primary Decomposition) We compute the primary decomposition of 

I={x^~y,y^~z,z^-x)eQ[x,y,z]. 



il 
12 
i3 
o3 



S = QQ[x,y, z] 

I = ideal (x' 2-y , y " 2-z " 2 , z " 2-x) 

dim I 





i4 : degree I 
o4 = 8 



15 : bpd = BPD I 
[ . . . ] 
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06 = { ideal ( z+ww_6-l , y-ww_6 + l , x+ww_6 ) , 

ideal ( z+ww_6, y+ww_6 , x-ww_6 + l ) , ideal (z + l,y-l,x-l), 
ideal (z-l,y-l,x-l), ideal ( z-ww_6 , y+ww_6, x-ww_6 + l ) , 
ideal ( z-ww_6+l , y-ww_6+l , x+ww_6 ), ideal(y,x,z"2)} 

17 : intersect bpd == sub (I, ring bpd#0) 

07 = true 

The function bpd is a shorthand for binomiaiPrimaryDecomposition, which can also be 
used in the long form and offers some options. The primary decomposition of / into bino- 
mial ideals exists in <[^{^(,)\x,y,z], so bpd created this cyclotomic field, calling the primitive 
sixth root of unity ww_6. Observe that the ideal has a double zero at the origin. In 17 we in- 
tersect the result to confirm that the decomposition is correct. The result of the intersection 
is defined over the extended polynomial ring ^(^e)!-^,^,^], and can be compared to / only 
after mapping it to that ring. 

This concludes our overview of the functionality of Binomials and we move on to the 
discussion of some large primary decompositions. 



2 A nonradical Commuting Birtli and Deatli Ideal 

In this section we study the commutative algebra of discrete time commuting birth and 
death ideals. One-dimensional birth and death processe s are among the simplest Markov 
chains that are considered in modeling random processes ( Latouche and Ramaswamill999l) . 



In the discrete time case, man y of their pr o pertie s can be derived from the explicit spectral 
theory of transition matrices. lEvans et al (1201 oh give motivation to consider generalized 



processes that correspond to Markov chains on multi-dimensional lattices, and as most of 
the one-dimensional theory does not apply there, the authors strive to identify subclasses 
with nice properties. The work suggests commuting birth and death processes which are de- 
fined by transition matrices having the property that transitions in the different dimensions 
commute. After reformulation, these conditions can be seen to result in binomial conditions 
on the entries of the transition matrices, that is, a binomial ideal. The toric component of 
this binomial ideal nicely relates to an underlying matroid as discussed in the paper De- 
termining primary decompositions of commuting birth and death ideals poses interesting 
challenges in combinatorial commutative algebra. 

Computational results given in this section tend to be very large. We have therefore 
stored them on a web page, which also contains additional scripts to reproduce the results: 

|http : / / personal- homepages . mis . mpg . de/kahle/cbd/l (15) 

We now define the binomial ideals under consideration. The ambient polynomial ring 
has variables corresponding to the edges of a regular grid. For fixed integers «i , . . . , n^, let 

m 

£:=n{0,...,«-l}, (16) 

be the usual ;7z-dimensional bounded regular grid with edges between vertices that differ by 
±1 in exactly one coordinate. Here it is sufficient to consider only the cases m = 2, 3. For 
each edge in the grid we define two variables, one for each direction. In the two-dimensional 
case the authors used the notation \[R,L,D,U] to denote a polynomial ring in the variables 

{Rij : < / < «i,0 < < 112] U {Uj : < ; < «i,0 < ; < in] U 
{Dij : < / < ni,0 < < 112} U {Uij : < i < riuO < j < 112} , 

where Rij is supposed to represent a right move starting at position ij and so on. In the 
case m = 3 one can, in a natural way, extend the set of variables by introducing letters 
F and B and three indices for each indeterminate. The set of commuting birth and death 
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processes is defined by the binomial equations (3.1) of lEvans et a] These equa 



tions arise in quadruples, coming from squares in the graph E, by which we mean in- 
duced subgraphs G of E that are isomorphic to the usual square. Denoting its vertices 
by {(m, v), (« + e/, v), (m, v + ey), (m + e,, V + ej)}, the corresponding ideal encodes that the 
two paths joining opposite vertices are equivalent: 

•= (U(u.y)R{u,v+ej) -R{u,v)U(u+ei,v)-, - R{uy+ej)D(u+ei,v+ej) j 

^{u+ei.y+ej)D(uy+ej) ^ D {^u+ejy+e j)l^(u+ej,v) i l^(u+eiy)U(uy) — L^(„+e,.,v)i'(i(+e,-,v+ej) ^ • 

(18) 

The commuting birth and death ideal is the sum of all l'^, where G runs through the induced 
squares of E. 

G square in E 



In the case m = 2,3 these ideals have been denoted /("i'"2) and /("i'"2'"3) bv Evans et al 
(l2010l) . 



Example 7 The graph E for m = 2 and n\ — 112 — I is just a square and /(' ') is generated 
by four binomials 

DiiLio -LiiOoi, L\qUoo-U\oL\i). 

If m = 3 and n\ = n2 = no, = \, E is the 3 -cube and the squares arise from facets. Thus, 
/( 1:1:1) js generated by 24 pure difference binomials, 4 for each facet. 

On the web page (flSl l one can download Python scripts that generate Macaulay 2 code for 
the rings and ideals in the cases m = 2,3. The following shows an example how to use the 
script imn .py on the command line to generate /(^'^): 

> . /Imn .py 2 2 

— Macaulay 2 Code for the Commuting Birth and Death Ideal: 

— m = 2 , n = 2 

S = QQ[ROO,UOO,R01,D01,U01,R02,D02,R10,L10,U10,R11,L11,D11,U11, 

R12, L12, D12, L2 0, U2 0, L21, D21, U21, L22, D22] ; 
I = ideal 

(UOO*R01-ROO*U10, R01*D11-D01*R00, D11*L10-L11*D01, L10*U00-U10*L11, 
U01*R02-R01*U11, R0 2*D12-D02*R01, D12*L11-L12*D0 2, L11*U01-U11*L12, 
U10*R11-R10*U2 0, R11*D21~D11*R10, D21*L2 0-L21*D11, L2 0*U10-U20*L21, 
U11*R12-R11*U21,R12*D22-D12*R11, D2 2*L21-L22*D12, L21*U11-U21*L22) ; 



In Evans et al the authors discuss the primary decompositions of /'^^'^\ /(IjIjI)^ and 



smaller examples. They state that these computations could not be carried out with the 
standard implementations, but were derived in an interactive session. The current imple- 
mentation of Binomials computed the 199 prime components of /(^'^^ in 100 seconds and 
took 123 seconds to decompose on the author's 1,6 GHz laptop. As mentioned be- 

fore, computing the minimal primes directly is even faster and can be completed in half of 
the time. 

Based on their results, Evans, Sturmfels, and Uhler conjectured 

Conjecture 1 For any grid E, the ideal is radical, its prime decomposition consists of 
pure toric ideals and is independent of the coefficient field. 

Here a pure toric ideal is an ideal generated by variables and pure difference binomials. 
Evans et a i (l201(ih 3rove that every associated prime of is a pure toric ideal. However, 



using Binomials we have derived the following counterexample to radicality. 
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codimension 


16 


17 


17 


18 


18 


19 


19 


20 


20 


21 


21 


22 


22 


# of components 


1 


14 


2 


107 


91 


356 


612 


527 


550 


212 


120 


38 


8 


gen. max degree 


1 


1 


4 


1 


6 


1 


5 


1 


4 


1 


2 


1 


3 


degree 


1 


1 


64 


1 


4012 


1 


144 


1 


36 


1 


12 


1 


3 


monomial 


y 


y 


n 


y 


n 


y 


n 


y 


n 


y 


n 


y 


n 



Table 1 Statistics on the primary components of sorted by codimension. Monomial components have been 
separated from binomial ones as indicated in the row "monomial". The row "gen. max degree" gives the maximal 
degree of a generator in this codimension while "degree" refers to the maximal degree among components. The 
toric component is generated in degree 6, of codimension 18 and degree 4012. 



It 


1 


2 


3 


4 


5 


6 


# of components 


3 


11 


40 


139 


466 


1528 



Table 2 Prime decompositions of/''"' 



Theorem 2 The ideal /P^^) is the intersection of 2638 primary binomial ideals whose prop- 
erties are given in Table\l] Among these are 10 components that are not prime, and thus 
/'^'^' is not radical. The 10 associated primes of these components are all embedded and 
of codimension 20. The radical Rad(/'-'^'^-') is the intersection of 2628 minimal primes and 
given by the following ideal: 

/(2,3) _|_ {DqiRot,RioLi2U2iL22D23 - C/oi^03i'10^13£'2li'23^23; 

U00R02R12L13L20D22U22 - RQQDo2RnLl3U2oU22L23, 

2 

RooUoiRo3Li(,Ri$U2oL23D23 - U01R03R13L13U20L23D23, 
ROOD02L10R13L13D21U22L23 - D02R03R13L13D21U22L23, 



(21) 



UQQR02R03Rul-\'il-2QD22D23 " ^Ooi'o2-Ro3-'^13^'13f^2oi'23£'23 , 
UQoRo3RloI^\2J^2(iU2\L22D23 — t^00^03^12^13t^21^22i'23£'237 
-'?00-D03-^'ll-'?13£^20^21-D22i'23 - UQ0R03L12R13L2QL22D22D23, 
RQ1U02L1QR11R13D21U21L23 -£>01^02^10^12il3t^21 1/22^23, 
D01R02R10R12L13L20D22U22 — DoiRo2Rl2Rl3Ll3D22U22L23: 
DmRQ3Rwl-\2l-\3U2\L2.2U22 - t^01^03i'10^13i'13£'2lf^22i'23)- 

One should note the two squares of variables in the third and fourth generator of Rad(/'^'^^). 
To produce these results one can use the functions bpd and binomiaiMinimaiPrimes. The 
author's computer determined the minimal primes in approximately 4 hours. Taking the 
intersection of these primes took another hour on a 2,8 GHz AMD Opteron. Care has to 
be taken when computing intersections of many primes. In Macaulay 2 versions 1.2 and 
below, using the command intersect directly on a large list of primes will not terminate. 
If one does the intersection manually with a loop, intersecting only two ideals at a time, 
everything is fine. Computing the cellular and primary decomposition was more delicate. 
It took several days and used about 5 GB of RAM. In fact, the original computation of the 
cellular decomposition was done with a slightly different algorithm which only works if the 
toric compo nent is isolated. We first computed the toric component T independently with 
the tool 4ti2 ( 4ti2 teaml2007 ) and then removed it by computing the saturation (/f^^^) : T"). 



The cellular decomposition of this ideal was easier to compute. Surprisingly this is not 
always the case. For some ideals /, with toric component T, the saturation (/ : r°°) is just 
too complicated to be computed with Macaulay 2. In some cases, simply doing the cellular 
decomposition with Algorithm[T]is faster 

To complete this computational study, we have also investigated the ideals for 
« < 6. It was not possible to find a counterexample there. 

Theorems The ideals n = 1,...,6 are radical. The respective numbers of prime 

components are given in Tabled 
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Concluding this section we find that the conjecture turned out to be false in full generality. 
It might however hold for the ideals and the associated primes could still be pure 

toric ideals for all . 



3 Conclusion and further directions 

We have presented algorithms for binomial ideals together with an implementation in 
Macaulay 2. It covers the case of pure difference binomial ideals, and it remains a future 
task to extend it to other cases, in particular to finite fields. 

A natural area for application of this software is the field of algebraic statistics, where 
analyzing the solutions of polynomial equations is of central importance. As mentioned 
in the introduction, describing conditional independence models is naturally connected to 
primary decomposition and also a very actively pursued research direction in algebraic 
statistics. The author hopes to facilitate experimentation with the availability of Binomials. 

Many operations that can be carried out with binomial ideals have been translated to 
operations on exponent vectors, or on the associated partial characters. By "making them 
combinatorial" significant speedups can be achieved. The computation of the associated 
primes is an example. Co mputing binomial primar y components is more delicate; the Hull 



operation is a bottleneck. Dickenstein et all (2008) give an explicit lattice point character 



ization of binomial primary components, but it seems not easy to use these results for 
computation. In the examples we have considered here, the Hull operation only marginally 
contributes to the total computation time. This is due to the fact that most of the compo- 
nents in our decompositions are prime ideals. In this case most of the computation time 
is spent on cellular decomposition, which in turn consists of many ideal saturations. Thus, 
from the author's point of view, software for binomial ideals would greatly benefit from a 
solution to the following problem; 

Problem 1 Develop a specialized algorithm to compute, for any (cellular) binomial ideal 
/, the "partially saturated" ideal 

^ . (22) 

The software 4ti2 implements the project-and-lift algorithm, a fast algorithm for computing 
the saturation 





(23) 

It seems natural to extended the program to solve the above problem, and Binomials is 
prepared to incorporate it upon availability. 

Finally, a natural approach to continue this work is to investigate decompositions that 
are finer than cellular decompositions, but not as fine as primary decompositions. In this 
direction one could aim at a separation of the combinatorial operations like cellular de- 
composition, and the field dependent operations like saturation of partial characters. The 
combinatorial operations should be connected to the combinatorics of the underlying semi- 
group ring. One can ask for the finest decomposition of a pure difference binomial ideal 
into pure difference binomial ideals, even if it is not primary. This might be interesting for 
applications where factorization of univariate polynomials is not of great importance. For 
example, if a component is generated by x^^ — I, we would like the algorithm to stop, since 
we know the result of this decomposition, and don't want the 19 cases to clutter up the 
output. It will be the subject of future work to investigate these possibilities. 
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